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The momentum spectra of particles in different high energy processes, such as e + e~ anni- 
hilation, pp and nuclear collisions in the limit p, pr — ► exhibit similar properties because 
of the dominant role of coherent soft gluon bremsstrahlung. We observe the following gen- 
eral features: the inclusive particle density approaches a limiting behaviour and becomes 
independent of primary collision energy; furthermore, it becomes proportional to the QCD 
colour factors Ca,Cf which appear in the Born term for the respective minimal partonic 
processes. In this limit, nuclear collisions reach with good accuracy participant ("wounded 
nucleon") scaling. Particle ratios in the low momentum region display a universal be- 
haviour. Future measurements at the LHC will provide crucial tests for the contributions 
from additional incoherent multi-component processes. 

1 Limiting soft particle emission 

We consider the invariant hadron density in the limit px — > or p — > and we define 

dN 



° d 3 p 



(1) 



If we calculate the momentum spectrum of gluons within perturbative QCD we find that the 
lowest order contribution, i.e. the Born term, dominates in the soft limit. This leads to some 
universal features for the various processes although they show very different features with 
increasing complexity in the multi-particle final states: quark and gluon jet production in 
e + e~, spectator jets and underlying event in pp and collective phenomena in AA collisions. 
Nevertheless, in the soft limit, we expect for the final state gluons the following properties: 

1. momentum spectra become energy independent; 

2. the relative normalisation of spectra in different processes is given by the colour factors which 
appear in the Born term for the respective minimal partonic process. 

The arguments hold for partons and we assume these properties are also valid for hadrons. 

First, we present qualitative arguments for this behaviour. A soft gluon is coherently emitted 
by all final state partons. Due to its large wavelength it cannot resolve the intrinsic structure of 
the bunches of partons and so it "sees" only the overall colour charge of the primary partons. 

* Speaker 



ISMD2010 



1 



The details of this study with the emphasis on pp and AA results are presented in Ref. [T] . 
The corresponding properties for e + e~ annihilation were obtained already some time ago [2] 
(for a recent overview, see [3]). 

The inclusive spectra can be derived using evolution equations, similar to the well known 
DGLAP equation, but with modification at small momenta to include coherence effects of soft 
gluons (for a review, see e.g. Ref. j^j). In the simplest double logarithmic approximation 
with fixed coupling a s one finds for the inclusive distribution in angle and energy k the 
perturbative expansion 

dN a 2C a 4Nc C a 2 E kr 
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The first term on the r.h.s. represents the well known Born expression for the soft brcmsstrahlung 
which is independent of jet energy E and holds the colour factor for the quark or gluon jet, 
respectively, with C a = Cf and C a = Ca- This term dominates if the transverse momentum 
kx ~ kO approaches the lowest value at the cut-off Qq, and this property remains valid in more 
exact calculations. The application to hadrons assumes similarity of parton and hadron spectra 
as required by the LPHD concept [3] with small Qo of few 100 MeV and accounting for the mass 
effects. Then a good description of the e + e~data in a wide range of energies s/s = 3 . . . 200 
GeV is obtained demonstrating the approach to a limiting behaviour for p — > 0. 

The difference between quark and gluon jets and the appearance of colour factors in the Born 
term ([2]) is best demonstrated in the measurement of soft radiation in 3-jet events perpendicular 
to the jet plane for which case the formulae for particle densities are given in [2]. Varying the 
inter-jet angles one can interpolate between the collinear configurations corresponding to qq 
and gg dipoles. This allows a test of the prediction 

I* 3 /If = C A /C F . (3) 

The measurement by DELPHI [6] confirmed this expectation and allowed a determination of 
the ratio in ^ as 2.211 ± 0.053 consistent with Ca/Cf = 9/4. In contrast, the ratio of total 
multiplicities in quark and gluon jets acquires large corrections to the ratio Ca/Cf- 



2 High-energy pp collisions 

We consider next the soft particle production in the non-diffractive "minimum bias" events in 
pp collisions. As in e + e~annihilation we look for the minimal partonic process which could 
be responsible for the very soft gluon brcmsstrahlung. We assume that the minimal process 
of lowest order corresponds to the semihard onc-gluon exchange between any two partons in 
the proton with a dominantly small scattering angle and a non-vanishing cross section at high 
energies. The exchange of a gluon between the partons in the proton leads to two separating 
outgoing partonic systems which give rise to gluon brcmsstrahlung from the effective colour 
octet dipole. In the simplest case each proton splits into a quark-diquark pair which scatter via 
one-gluon exchange. Also more complex partonic processes could occur with two colour octet 
systems, as discussed in Ref. pQ. 

Therefore, we expect limiting behaviour of the soft particle density in pp collisions, /q P , for 
p — > or pt — > at fixed rapidity y. Furthermore, we predict the ratio to the corresponding 
density in e + e~ annihilations 

p^O: IF/I?*' hC a /C f , (4) 
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Figure 1: Soft limit I pp of the invariant density Ej^ of charged particles [{h + + h )/2] in pp 
collisions as a function of c.m.s. energy i/s (for exponential extrapolation, normalised to a in ). 

similar to the case of primary gg and qq dipolcs. The px distributions in both processes differ 
and, therefore, the integrated multiplicities differ as well. 

In order to test these predictions we studied the invariant cross sections E-^ measured in 
the energy range y/s = 20 . . . 1800 GeV at the colliders at BNL, CERN and Fermilab. After 
normalisation to the inelastic cross section, ai n , we found the corresponding density Iq from 
extrapolation of the pt spectra to pr = 0. Most experimental groups fitted their data using 
parametrisations with the small pr behaviour 



which work well down to the smallest measured values of pr ~ 0.1 GeV. We derive Iq P = A/ai n 
from these fits. The functional form ([5]), however, is not analytic at px = 0. There is no such 
problem from the theoretical point of view for the "thermal" parametrisation in terms of the 
transverse mass mr instead of px 



This form was applied by the PHOBOS collaboration [7], and they obtained a good fit to their 
AuAu data down to px ~ 0.03 GeV. With such extrapolation the resulting value of Iq P is found 
~25% lower as compared to the fit with parametrisation ([5]). 

The results from the available published exponential extrapolations and normalisation to 
inelastic cross sections Ui n are shown in Fig. [TJ The data on I pp are seen to be consistent with 
energy independence over two orders of magnitude in energy with fluctuations around a 
mean value Iq P sa (7 ± 1) GeV -2 . The different functional forms of extrapolations mentioned 



E-pr- = Acxp(Bp T + ...), 



(5) 
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above do not affect this conclusion pQ. It should be noted that the pT-intcgrated rapidity 
density ^ doubles in this energy range. 

The observation of energy independence is similar to that in e + e~annihilation and corre- 
sponds to the expectation from coherent bremsstrahlung. In contrast, for an incoherent super- 
position of several sources one would expect a rising Iq P , since the number of such processes 
would typically increase with energy It will be interesting to obtain corresponding results from 
the LHC at the higher energies to test this possibility. Some results were obtained by CMS [8] 
which show the convergence of spectra at different energies towards pr = 0. 

Since the soft gluon emission is driven by the lowest order diagram we predict the same 
limiting behaviour (da / dycPpT const in yfs at pr — > 0) to be observed at any rapidity and 
not only at y = and we expect a flat rapidity plateau in this limit. This is true, as long as 
we can neglect the kinematical boundary effects for large rapidities y. To study the rapidity 
dependence it is important to perform calculations with the exact kinematics using particle 
masses. The use of pseudorapidity r] = — lntan(#/2) will lead to a different spectral shape in n 
and pt as discussed in detail in Ref . [§] . 

For comparison of particle densities in pp and e + e~collisions in order to test eq. (0]) it is 
better to use for normalisation the non-diffractive cross section which we take as being 15% 
lower than the inelastic one, and we obtain lQ P nd ~ (8 ± 1) GeV _2 [T]. Then we find for the 
non-diffractive (minimum bias) events for a thermal fit the reduced value 

7 pp « (6±1) GeV" 2 . (7) 

In order to compare with e + e _ annihilation let us consider two results. 

1. The TPC/27 collaboration [10] directly compared their own data on e + e~annihilation 
with the pp data by the BS collaboration |llj . The pr spectrum in the pp collision is found to 
be steeper than that in e + e~ annihilation , which is determined with respect to the sphericity 
axis. The extrapolation down to the smallest pr yields a larger density in pp collisions by 
factors 2.0-2.7 depending on the type of the fit. 

2. Other e + e~ experiments presented cross sections as a function of particle energy. Using 
their exponential fits in this variable in the energy range y/s = 10 ... 29 GeV we find Iq e sa 
(3.3 ± 0.5) GeV~ 2 or the ratio 

I^/lf e ~ « (1.8 ± 0.4) ^ (2.4 ± 0.5), (8) 

where the first (preferred) number refers to the thermal and the second to the exponential 
extrapolation. This result is consistent with our expectation for this ratio Ca/Cf — 2.25. The 
approximate energy independence of the quantity I in both collision processes and that their 
ratio is close to 2 are remarkable and represent a serious argument in favour of the dominance 
of the elementary bremsstrahlung process, also in soft pp collisions. 



3 Low pr spectra in nuclear collisions 

For nuclear collisions the data are often presented in terms of a corresponding ratio to the pp 
collision according to the two limiting cases: 

1. We can define the "nuclear modification factor" 

n N coll 1 dN AA /dp T 

N coll dN pp /d PT ' {J> 
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which is equal to unity in the case of a pointlike interaction. Here the number of nucleon-nucleon 
collisions N co tt can be calculated from the Glauber model. 

2. Another limiting case applies to the bulk particle production without any hard scale. 
In such a case we expect that soft particles are produced coherently over a certain range r 
characteristic for a nucleon (or even nucleus) which results in a reduced rate. Indeed, all 
collaborations at RHIC [7J Q21 EH E] find that the ratio fi^j" is falling below unity for small 
Pt- Data in this kinematic region are conveniently presented by normalising to the number of 
"participating nucleons" 

R N part 1 dN AA /dp T 

(AW/2) dN pp /d PT ■ [ Uj 

Bialas, Blcszynski and Czyz |15j have introduced the concept of participating nucleons, or 
"wounded nucleons" in their terms, and observed that R AA rt ~ 1 for soft production, which 
means that each interacting nucleon should be counted only once and rescatterings be disre- 
garded. 

As before, we consider first the energy dependence of the particle density at px — > 0. The 
variation of px spectra between the two energies 62.4 and 200 GeV was studied for different 
centralities by PHOBOS [7J. The nuclear modification factors R AA " approach about the same 

values for the two energies at small px ~ 0.2 GeV, and the same is true for R A J A art . This implies 
that the energy dependence of AA data is as weak as that of pp data. 

Next we investigate the normalisation of px spectra in AA collisions. This follows from 
the observation [7] R A f A aTt — > 1 for small px, i.e. approximate validity of "participant scaling". 
We studied the issue further by combining the high accuracy pp and AuAu data from STAR 
[12"] at 200 GeV with the corresponding nuclear data from PHOBOS that extend towards low 
Pt ~ 0.03 GeV. Using the thermal paramctrisation ^ we find 

I^/If* « 160 ± 17, (11) 

which is compatible with the calculated N part /2 = 172 (±15%). Therefore, the normalisation 
in the soft limit is consistent with 

PT^O: R^F-H and I AuAu « ^p± i». (12) 

We find it quite remarkable that this simple relation works with good accuracy of about 10% 
which corresponds to the accuracy of both the measurement and the theoretical calculation. 
Note, in particular, that relation (|12j) is only valid in the limit pt — > and is violated by about 
50% already at px ~ 0.5 GeV (see STAR [M]); also note that the observed particle density for 
the "wounded nucleon" interactions is strongly suppressed (ex Np £ Tt ) by the coherence effects 
and it is about six times smaller than for the incoherent superposition of nuclear scatterings 
(oc Neon). 

Within our approach we should relate this normalisation of particle density in the soft limit 
(|12l) to the minimal parton configuration responsible for the emission of a large wavelength 
gluon. In this limit with "participant scaling," nucleons of one nucleus interact with some 
nucleons of the other nucleus several times but the outgoing partonic system acts again like a 
colour octet source such as if it had scattered only once as in pp collisions. One could think of 
two mechanisms leading to such a result. 

The primary interaction in a single nucleon-nucleon collision responsible for a "minimum 
bias" event leads to the break up of the scattered nucleons dominantly into quark-diquark 
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systems which form colour singlet or octet states only. One could argue that such a semihard 
collision with px ^5 1 GeV may not resolve higher Fock states in the proton of typically smaller 
size. Then the soft hadrons will not be sensitive to higher partonic excitations either^ Another 
argument can be based on an assumption that at high energies the multi-gluon exchange related 
to multiple scatterings of the same nucleon dominantly proceeds via the colour octet exchange, 
while the higher colour multiplets (such as decuplet or 27-plet) are disfavoured. Note that in 
the non-Abelian gauge theory the intercept of the reggeised gluon ac(0) = 1 (i.e. the spin 
of the on-mass-shell gluon is equal to one) while the j-planc singularities corresponding to the 
exchange of the higher colour multiplets are situated to the left of j = 1 [TH]. Therefore, the 
exchange of higher colour multiplets are asymptotically suppressed. It remains a challenge to 
understand the rapid disappearance with pt of the limiting density f| 12[) . 

Another interesting simplified structure appears for particle ratios in the low pt region. 
If we compare [T] the pt dependence of such ratios from e + e~annihilation (j>t with respect 
to sphericity axis) with those from pp scattering [TU] and also ratios from pp and AA colli- 
sions by PHENIX and STAR we observe the convergence of these spectra in the region of low 
Pt 0.5 GeV while they are quite different otherwise. This could be related to the universal 
bremsstrahlung from the colour sources in the considered processes in this soft limit. This 
observation implies that the soft particles decouple from equilibration in nuclear collisions. The 
soft particles in the central rapidity region are produced first and they stay behind if the system 
expands. 
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A model based on wounded quarks and diquarks has been developed in Ref. |17j . but for the description of 
p^-integrated rapidity distributions. 
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